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Abstract 

We use gauge theory/string theory correspondence to study finite temperature critical 
behavior of mass deformed M = 4 SU(N) supersymmetric Yang-Mills theory at strong 
coupling, also known as M = 2* gauge theory. For certain range of the mass parameters, 
Af = 2* plasma undergoes a second-order phase transition. We compute all the static 
critical exponents of the model and demonstrate that the transition is of the mean- 
field theory type. We show that the dynamical critical exponent of the model is z = 0, 
with multiple hydrodynamic relaxation rates at criticality. We point out that the 
dynamical critical phenomena in M = 2* plasma is outside the dynamical universality 
classes established by Hohenberg and Halperin. 
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1 Introduction 

Gauge theory/string theory correspondence [1] presents a solvable framework to study a 
large class of strongly interacting four-dimensional gauge theory plasmas. In a nutshell, 
the solvability of these models comes from ability to approximate a dual string theory 
with a corresponding classical supergravity. Unfortunately, real QCD is not any one of 
the models studied. It is possible to reach QCD as a particular limit in some of these 
models, but in doing so the truncation of the string theory to a supergravity sector 
becomes inconsistent. Instead one attempts to discover common/universal features of 
strongly coupled gauge theory plasmas, and hopes that real QCD is in the universality 
class of the models studied. Typical examples of such universal properties are the 
strongly coupled plasma shear and bulk viscosities: 
■ the shear viscosity r\ to entropy density s ratio [5-8] 



■ the bulk viscosity bound [9] 

$*«(H). 

Although above properties of strongly coupled plasmas have been observed (or in case 
of the shear viscosity derived) in holographic setting, it is not clear why and how 
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these universalities arise, or how to properly define the corresponding universality 
class: while the shear viscosity ratio in universal in 2-derivative supergravity 1 (or a 
phenomenological model of thereof), it can be violated in full string theory [10-12]; 
while the bulk viscosity bound is satisfied in all models of supergravity derived from 
string theory, it can be violated in some phenomenological models of gauge/gravity 
correspondence [13]. 

A more common notion of the 'universality' arises in the theory of continuous crit- 
ical phenomena. In this paper we follow up the work of [14-16] and focus on static 
and dynamic properties of strongly coupled non-conformal gauge theory plasma in the 
vicinity of the second-order phase transition. In [16] it was noticed that there was a 
tension between the hyperscaling relation among the static critical exponents at the 
second-order phase transition, and the expectation that in the planar limit the transi- 
tion should be of the mean-field type, i.e., with vanishing anomalous critical exponent 2 . 
Direct computation of critical exponents for the second-order phase transition in M = 4 
SYM plasma at finite temperature and the chemical potential for a global U(1)r charge 
confirmed the vanishing of the anomalous critical exponent. Further, the dynamical 
critical exponent of this transition was shown to be z = 4, even though the background 
geometry at criticality did not exhibit az = 4 Lifshitz-like scaling. In other words, the 
transition detailed in [16] explicitly showed that the dynamical scaling properties can 
be "emergent" and should not be necessarily "enforced" on the background geometry 
of the holographic dual. 

Although we restrict our attention here to a second-order phase transition in mass- 
deformed Af = 4 SYM (also known as Af = 2* gauge theory [2-4]), we emphasize that 
the holographic (static) universality class of this transition includes also a cascading 
gauge theory [22]. In section 2 we review the holographic duality for Af = 2* gauge 
theory plasma. Critical phenomena in Af = 2* plasma from both the gauge theory and 
the dual gravitational perspective is discussed in section 3. Some of the static critical 
exponents of the second-order phase transition in this plasma, namely {a, /3, 7, 5}, were 
computed in [15]. We directly compute the remaining static critical exponents {u, r/} 
and the dynamical critical exponent z of the theory in section 4. We collect all the 
results in section 5. 

lr rhis translates into an infinite t' Hooft coupling limit on the gauge theory side. 

2 Contrary to some statements in recent literature (as in [21] for instance), we take a perspective 
here that for a second-order phase transition to be of a mean-field type the anomalous critical exponent 
must vanish — whether or not the other critical exponents are integers or not is irrelevant. 
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2 J\f = 2*/PW holographic duality 



In this section we briefly review the main features of the holographic duality between 
Af = 2* SU (N) gauge theory and the Pilch- Warner (PW) geometry of type IIB super- 
gravity. We refer the reader to the original work for further details [2-4, 17-20]. 

Consider maximally supersymmetric Af = 4 SU (N) Yang-Mills theory in the planar 
limit (g Y M ~ * 0) N —> oo with A = g Y M^ kept fixed) and for large 't Hooft coupling 
A ^> 1. According to Maldacena correspondence [1] this superconformal theory is 
equivalent to a classical type IIB supergravity on AdS 5 x S 5 . A duality between a 
SYM and a supergravity can be extended (on both sides) away from the conformal 
point [2-4]. On the gauge theory side, a massive deformation of Af = 4 superpotential 



WV=4 = ^- Tr(b,gU) , (2.1) 
9ym vl J 7 

where {Q, Q, $} are Af = 1 adjoint chiral superfields, to 

W M=4 -)■ WV=2* = W M=4 + — — (TtQ 2 + TtQ 2 ) , (2.2) 

9ym v j 

breaks half of the supersymmetries. This mass-deformed theory is known as Af = 2* 
gauge theory. When m/0, the mass deformation lifts the {Q, Q} Af = 2 hypermul- 
tiplet moduli directions, resulting in (N — 1) complex dimensional Coulomb branch 
parametrized by 

N 

$ = diag(ai,a 2 , • • • ,a N ) , y~] a, = . (2.3) 

i=i 

We study Af = 2* gauge theory at a particular point on the Coulomb branch moduli 
space [3]: 

OiG [-ao,ao], al = m2g2YMN , (2.4) 

7T 

with the (continuous in the large A^-limit) linear number density 

2 / f a " 

p(a) = ^ r z-Ja 2 -a*, / da p(a) = N . (2.5) 

m 9 Y M V J-ao 

The reason for such an esoteric choice for a vacuum of the theory is simply because 
we know a dual holographic description of the theory (as a Pilch- Warner geometry [2]) 
only at this point [3]. Extending the correspondence to the rest of the moduli space is 
an important unsolved problem. 
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Notice that the deformation (2.2) is actually a deformation of a CFT by two different 
operators: a dimension-2 operator (a mass term for the bosonic components of the 
{Q, Q} hypermultiplet) and a dimension-3 operator (a mass term for the fermionic 
components of the {Q,Q} hypermultiplet). According to AdS/CFT dictionary [23], 
a scalar gauge-invariant operator of dimension A is dual to a scalar field of mass 
m 2 L 2 = A(A — 4) of the five- dimensional dual gravitational description. These two 
mass- deformation operators are the a and x scalars of the Pilch- Warner effective action 
[2]: 

S = f deV^ £5 = -r^r [ d?y/=g [\R - 3(da) 2 - (d X ) 2 - V] , (2.6) 
where the potential 3 



1 

V = — 
16 



1 fdWY fdW 



3 \ da J \ dx 
is a function of a and Xi an d i s determined by the superpotential 



lw 2 , (2.7) 



W = -e~ 2a - -e 4a cosh(2x) . (2.8) 
2 

In our conventions, the five- dimensional Newton's constant is 

In what follows we focus on equilibrium thermal states of M = 2* plasma. Their 
holographic dual is represented by a regular black brane solution in the effective action 
(2.6) [17,19] 

dsl = e 2A (-(1 - xfdt 2 + dx\ + dx\ + dx 2 3 ) + g xx dx 2 , (2.10) 

with g xx = g xx (x), A = A(x), a = a(x) and x — x( x ) being functions of the radial 
coordinate x G [0, 1] only. Note that x — > + corresponds to the asymptotic AdS 5 
boundary, while x — > 1_ to a regular Schwarzschild horizon. The temperature and 
the mass parameters of the plasma are encoded in the asymptotic behavior of the 
supergravity fields {A,a,x}- Specifically, near the AdS$ boundary we have 

e a = p = 1 + x 1 / 2 (p 10 + Pll \nx) + --- + x k / 2 Pki WxJ +••• , (2-H) 



3 We set the fivc-dimensional gauged supergravity coupling to one. This corresponds to setting the 
radius L of the five-dimensional sphere in the undeformed metric to 2. 
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X = Xox 



1/4 



1 + x 1/2 (Xio + Xu In a;) + • • • + a; fe/2 ^ X* + 
1/2 (aio + an lnx) + • • • + x k/2 I a ki W x ) 



a = x 



and 



+ p fe (l-:r) 2fe + - 
. + Xk (l-x) 2k + 



P = ph + Pi(l - a;) 2 + • 

X = Xfe + xi(i - + 
a = a h + ai(l — x) 2 H h a fc (l - x) 2k + 

near the regular Schwarzschild horizon. In (2.13), (2.16) a(x) is defined as 

A(x) = In <5 3 — - ln(2x — x 2 ) + a(x) . 



(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 

(2.17) 



In was shown in [19] that given {<5 3 , pn,Xo}, there is a unique singularity-free solution 
of (2.6) representing the equilibrium state of M = 2* plasma. On the gravity side, the 
coefficients of the leading asymptotics, namely {5$, pu,Xo}, determine the remaining 
6 parameters (2 in the UV and 4 in the IR) of the solution: 

UV : {pio, Xio} , 

IR : {Ph, Xh, ah, a>i} . 
It is possible to unambiguously relate the gravitational and the gauge theory data [19]: 
■ the plasma temperature T, and the masses {m^, rrif} of the bosonic and the fermionic 
components of the M = 2 hypermultiplet are given by 



(2.18) 



2tt 



P = U^pne^ , 'f = 2 3 /V X(1( 



3a h 



(2.19) 



■ and the plasma free energy density J 7 , the energy density £, and the entropy density 
s are given by 



si 



32nG, 



1 + p 2 n (24 - 96 In 5 3 + 24 In 2) - 2Ap w p u + 2 X 2 oXlo 



4 2 



+ Xol 9 " 3 ln2+ 3 ln<5 3 



(2.20) 



£ = T 



1 ~ ffip'ia.h 
J—S 4 S ~ -2- 

krG 8 3 ' 4G 5 



To recover the Af = 2 supersymmetric PW vacuum (2.4), (2.5) we need to set T = 
and fine-tune the masses m h = rrif = m. 
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Figure 1: (Colour online) The dimensionless temperature (left plot) and the speed 
of sound Cg (right plot) of the strongly coupled M = 2* plasma with mf = and 
m& 7^ as a function of the dual gravitation parameter p n . 

3 Critical phenomena in M = 2* plasma 

The phase diagram of M = 2* plasma was studied in details in [19,20]. It was found 
there that whenever mj < ml, the theory undergoes a second-order phase transition, 
with the critical temperature T c = T c (j^j • All these transitions are in the same 
universality class, and thus we can restrict our attention to mj = 0, m& ^ case: 

m f = 0: ^» 2.32591. (3.1) 

-L C 

We now recall the main characteristics of this transition [15]: 

■ The left plot on Figure 1 represents the dependence of the dimensionless temperature 
on the gravitational parameter pu, see (2.19). The transition is associated with the 
minimal accessible temperature, to be identified with T c , in the plasma for isotropic 
and homogeneous equilibrium state 4 . For each temperature T > T c there are two 
phases — the "ordered" phase (blue curves), and the "disordered" phase (red curves). 
The right plot on Figure 1 represents the square of the speed of sound cj. as a function 
of pu. Notice that the hydrodynamic modes in the "disordered phase" are unstable 
(as c 2 s < 0), and thus must condense. It is tempting to conjecture that the equilibrium 
state in the plasma at T < T c breaks translational invariance, and represents the end 
point of condensation of hydrodynamic modes [24] . 

4 This feature of the transition is also observed for the phase transition in Af = 4 SYM plasma with 
a single U(l) C SU(4)r R-symmetry chemical potential [16]. 
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Figure 2: (Colour online) Free energy densities fl Q of the "ordered" phase (blue curves) 
and fid of the "disordered" phase (red curves) as a function of pn (left plot) and 
(right plot) of M = 2* plasma with rrif = 0. 

■ The free energy densities of the stable "ordered" phase (blue curves) Q D and the 
unstable "disordered" phase (red curves) Qd as a function of the gravitational parameter 
pn (left plot) and the dimensionless temperature ^ (right plot) are represented in 
Figure 2. 

■ It is convenient to recast the critical behavior in M = 2* plasma in that of a 3- 
dimensional ferromagnet. The thermodynamics of the latter is described by the Gibbs 
free energy W = W(t, "H) which depends on the reduced temperature t = (T — T c )/T c 
and the external magnetic field H. Once we identify 

w = n -n d , n = m b , (3.2) 

and introduce 

Ap u = p n -p c n , |A Pll | cxt 1 / 2 , p c n = 0.035187(6) , (3.3) 
we can compute the standard static critical exponents {a, (3,^,5}: 
'd 2 W\ ' blue 



-a_ S 



oc c 2 



red 



blue 



oc (Apn)- 1 oc r 1 ' 2 



red 



1 

a= 2 



oc -|Apu| oc -t 1/2 (3 = -, 



(3.4) 
(3.5) 



(&£f) oc |tr oc -^M oc d// 2 oc r 1 / 2 



\r = ( ^777 ) x IT *• -OiM -x (J,!* 1 x r J 2 =J> 7 = - , (3.6) 



8 



M(t = 0) oc \U-U C \ 1/S oct 1 / 5 oct 1/2 5 = 2. (3.7) 

Thus, 

{«,/3,7,<5}=(i 1, 1, 2) . (3.8) 



v 2 2'2' 

The remaining two critical exponents {u, 77} are more difficult to extract as they are 
related to the scaling properties of the magnetization two-point correlation function at 
criticality: 

G (^(W(5)>c^|^, (3,) 

{ e -\m 1 4 
with fccltl-", (3.10) 
| r1 -3 + 2-^ f = Q 

where £ is the correlation length. Under the static scaling hypothesis, 

W(t, U) = A~ 3 W(\ yT t, \ yn U) , G(q, t, H) = \ 2y ^ 3 G(\q, \ yT t, \ yH V) , (3.11) 

where i/t and y-^ are the two independent critical exponents; G is a spatial Fourier 
transform of (3.9). The static scaling hypothesis implies 4 scaling relations between 
{a, (3, 7, 5, v, 77}. In particular, using two of these relations 

2-a = 3u, 7 = 1/(2-77), (3.12) 

and (3.8), we find 



1,1^. (3.13) 



scaling 

Much like in [16], the non- vanishing of the anomalous critical exponent 77 conflicts with 
the expectation that for large-N gauge theory plasmas the continuous phase transitions 
in holographic models should be of mean-field type ri mean _f ield = 0. 

A relaxation of the system to equilibrium in the vicinity of the critical point is 
commonly discussed within the theory of the dynamical critical phenomena developed 
by Hohenberg and Halperin [25]. According to [25] a model is designated to a specific 
universality class based on the dimensionality, symmetries of the order parameter, the 
presence of any conserved densities, and any other properties that affect the static 
critical behavior. A representative of a given dynamical universality class is then 
characterized by a dynamical critical exponent z. This critical exponent determines 
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the scaling of the non-equilibrium (time-dependent) two-point correlation function of 
the order parameter (magnetization in our case) at criticality, i.e., 

G(u, q, t, U) = \ 2y ^- 3+z G(\ z to, Xq, X VT t, \ m H) , (3.14) 

for its space-time Fourier transform. The equilibration of a dynamical system is thus 
characterized by a relaxation time r 

rocT, (3.15) 

which (for z ^ 0) diverges at criticality. The absence of any conserved order parameters 
puts M = 2* plasma in the universality class of 'model A' according to the classification 
of Hohenberg and Halperin, and predicts 

= 2 + cr], (3.16) 

prediction 

where the constant c can be computed via renormalization group calculations in p = 
4-e,e < 1, spatial dimensions, and r\ is the anomalous critical exponent. 

In the rest of this section we introduce dynamical susceptibility of M = 2* plasma 
and explain how is can be used to compute the static critical exponents {is, r/} and the 
dynamical critical exponent z of its phase transition. 

3.1 Dynamical susceptibility of J\f = 2* plasma — gauge theory perspective 

Both the critical exponents {z/, 77} and the dynamical exponent z can be extracted from 
the dynamical susceptibility of the system. Consider the response of the system to the 
time-dependent inhomogeneous variations of the external magnetic field "H, 

d 3 k f dtu 



H^H + 5H(t, x), 5U = j j^- 3 j ==- e^~^ % uS . (3.17) 

At the linearized level the variation of the external magnetic field would produce a 
corresponding variation in the magnetization 5A4(t,x) ( M. w g for the Fourier compo- 
nents). Following [25] we introduce the dynamical susceptibility as 



M. 



, — . X^k = XT = 



dM 



(3.18) 



By the equipartition theorem, the static susceptibility 

Xfe = X w= o,fe> ( 3 - 19 ) 
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is related to the Fourier transform G(k) of the magnetization variation two-point cor- 
relation function 

G(x) = (5M(x)5M(d))sn=o , (3.20) 

as 

G(k)=T Xi . (3.21) 

Given the near-critical behavior of the correlation function (3.20) (see (3.10)), (3.21) 
implies that the static susceptibility x% nas a P°l e a ^ 

k 2 oc -r 2 , (3.22) 

in the vicinity, but not right at the critical point. On the other hand, right at the 
critical point 

Xfe -oc|£|- 2+ ". (3.23) 

The theory of dynamical critical phenomena [25] predicts that in the vicinity of the 
continuous phase transition, and for \k\ ~ the full dynamical susceptibility x^k 
will develop a pole at 

u ~ -if* , (3.24) 

with z being the dynamical critical exponent of the system. The frequency in (3.24) 
(in the hydrodynamic limit) defines a relaxation time r as 

t' 1 = iu <x £~ x . (3.25) 

To summarize, following the position of the poles in the static susceptibility Xk as 
a function of the reduced temperature t ^ 



1 

= — 



fc # 2 (*)cx-r 2 cx-t 2v , (3.26) 

\k\ 2 =k 2 ,{t) 



would determine the critical exponent v\ the critical exponent r\ is determined from the 
static susceptibility scaling at critical temperature, i.e., t = 0, as in (3.23). Likewise, 
scaling of the pole in the dynamical susceptibility x^ % * n ^ ne hydrodynamic limit as a 
function of the reduced temperature t ^ determines the dynamical critical exponent 

z: 

= ^ 



X u ,k 



=> iu*(t) oc t~ xv . (3.27) 

(u>=w»(t),fc-»0) 
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3.2 Dynamical susceptibility of M = 2* plasma — gravity perspective 

In case of Af = 2* plasma we identify the external magnetic field H with the bosonic 
mass m b , (3.2). The variation "H w £ would correspond to the variation in 

m b ->m b + 6m b (uj, k)e il£ ~ iuJt , (3.28) 

which on the gravity side can be induced by the variation in the non-normalizable 
coefficient pn the supergravity scalar p: 

Pn ->■ Pn + 5pu(w, k)e ik ' £ ~ iut , 5p u (w, k) oc 5m b (oo, k) . (3.29) 

The variation 5pu(uj, k) would produce a linearized response in the normalizable coef- 
ficient pio(u>,k) of the supergravity scalar p: 

Pn -+ Pn + S Pii (uj, fye**-™ p w -+ p w + S Pio (uj, k)e^ s -^ . (3.30) 

— # 

Thus it is natural to identify the variation 5pio(u>, k) with the variation in the magne- 
tization M.^ g 

5p 10 (u,k)(xM U! j;. (3.31) 
Finally, the dynamical susceptibility (3.18) is related to the dual gravitational data as 

8p 10 (u,k) . . 

X u ,k oc t • (3.32) 

opu{u}, k) 

The identification (3.32) is equivalent to the one made in recent analysis of the holo- 
graphic critical phenomena [14,16]. 

The holographic computation of the susceptibility (3.32) necessitates the analysis 
of the linearized fluctuations in the gravitational background (2.10). The relevant 
fluctuations were studied previously in [20,26]. We briefly review the basic setup of 
such computations here. Without the loss of generality we can assume that 

fc* = qS i 3 . (3.33) 

The linearized on-shell fluctuation of the gravitational scalar 5 a = In p, 

a^a + 4>{x) e iqX3 - iut , (3.34) 



5 For a critical phenomena with m,f =0we can consistently truncate the effective action (2.6) to 

x = o. 
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couples to the on-shell fluctuations in the background metric 

V^V + VW e ^-™\ (3.35) 

which form a helicity-0 representation with respect to rotations about rr 3 -axis: 

\htt) ht X3 , h aa h XlXl + h X2X2 , h X3X3 } . (3.36) 

Note that we partially fixed the background diffeomorphisms with 

htx = h X3X = h xx = . (3.37) 

Following [26] we introduce the diffeomorphism-invariant linear combinations of fluc- 
tuations 

Z H = A^-H tz + 2H ZZ - H aa (l + + 2 4( X - x )* H tt > ( 3 - 38 ) 

2{\ng XlX J 

where gu(x) and g xixi (x) are the corresponding components of the background metric 
(2.10), the derivatives are with respect to x and 

htt 9ttHtt j ht z gxix\Ht z , h aa gx\x\H aa , h X3X 3 9xixiH X3 x3 ■ 

(3.40) 

Introduce 

The equations of motion for {Z H , Z^} take the form 



— Z'h + C\\ Z' H + C12 ^ + C13 Z u + C14 , 
=Z'(j, + C21 Z# + C 2 2 z'^ + C 23 Z H + C 2 4 , 



(3.42) 



where the coefficients Qj are nonlinear functionals of the background fields {p, a} with 
explicit dependence on x and {tt>, q} [26]: 



C =C 



{p,a}; {m,q} 



(3.43) 



Since the equations (3.42) are homogeneous, we can always set the non-normalizable 
component of Z^ — which is the diffeomorphism-invariance analog of 5pn(cj, k) — to 
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Figure 3: (Colour online) The scaling (blue dots) of the inverse of the static suscepti- 
bility Xm=o,q=o in the vicinity of the critical point. The solid red line is a quadratic fit 
to data, the dashed green line represents p n = p\ x . 

one; the dynamical susceptibility is then proportional to the normalizable component 
of Z^. 

We can summarize now the boundary value problem whose solution would deter- 
mine the dynamical susceptibility Introducing 6 



z H = O(x) , = (In x + Z(tt>, q))x 1/2 + 0(x\n 2 x) , as x -)■ + . 

The normalizable component Z of z$ near the boundary is proportional to the dynam- 
ical susceptibility: 



4 Critical exponents {u, 77} and z of J\f = 2* plasma 

In previous section we explained how the dynamical susceptibility Xro,q can be used to 
extract the static critical exponents {^,77} an d the dynamical critical exponent z of 
6 The ro— and q— dependent rescaling are for convenience in further analysis. 




Xm, q oc Z(n>,q) . 



(3.46) 
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Figure 4: (Colour online) Poles of the static susceptibility in the vicinity of the critical 
point: Xwio,q=q* = 0- The solid red line is a quadratic fit to data, the dashed green 
line represents p n = p± v 

a phase transition. We also related this susceptibility to the normalizable component 
Z(yo, q) of the gravitational scalar, which non-normalizable component played the role 
of the (time-dependent and inhomogeneous) variation of the external magnetic field, 
see (3.46). Here, without going into the technical details of the analysis of the boundary 
value problem (3.42)-(3.45), we present the results. 

Figure 3 shows the inverse of the static susceptibility at q = (blue dots) in the 
vicinity of the critical point. The solid red line Zj it = Zj it (pn) is the best quadratic 
fit to the data. The vertical green line denotes the critical value of pn, i.e., p\ x (3.3), 
corresponding to critical temperature T c , see (3.1). The red line intersects the pn axis 
at p* 1; such that 



Z~ x 



= 

Pll=p*l 



Pn_ 1 



9.(3) xl(T 6 , (4.1) 



in excellent agreement with the critical behavior of xt deduced from the thermody- 
namics (3.6): 

Zjl oc Apu oc t 1/2 xt 1 °c t 1/2 . (4.2) 

Figure 4 presents the poles (blue dots) of the static susceptibility at q = q* in the 
vicinity of the critical point: 

X ro io,q=q, = . (4.3) 
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Figure 5: (Colour online) The scaling (blue dots) of the inverse of the static suscepti- 
bility Xre=o,q & t the critical point, p n = p\ v The solid red line is a quadratic fit to the 
data. 

The solid red line represents the best quadratic fit to the data, and the vertical green 
line denotes the critical value of pn, i.e., p xl , see (3.3). Notice that in the stable 
phase, i.e., for p n < p xl , in the vicinity of the phase transition the poles in the static 
susceptibility are for purely imaginary momenta, which implies the exponential decay 
of the magnetization density two-point correlation function (3.20). Furthermore, from 
(3.26) we identify the correlation length as 

(2ttT c 2 oc -q; 2 oc — 1—r oc +t~ 1/2 , < p\ x - p n < p c u , (4.4) 

|Apu| 

where we used the results of the fit and the relation between p u and the reduced 
temperature t (3.3). From (4.4) we extract the (static) critical exponent v: 

i oc r u oc r 1/4 v = - . (4.5) 

Given that the static critical exponent a = |, (4.5) implies that the hyperscaling 
relation (3.12) is violated 

2- a ^3v. (4.6) 

Figure 5 shows the inverse of the static susceptibility as a function of q (blue dots) 
right at the critical point p X \ = p\ x . The solid red line represents the best quadratic fit 
to the data 

(Zf*)- 1 = -5.55084 • 10" 6 - 0.30963 q 2 - 0.28859 q 4 + C(q 6 ) . (4.7) 
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Figure 6: (Colour online) The inverse of the dynamical susceptibility of = 2* plasma 
in a stable phase (blue curve) and an unstable phase (red curve) at q = 10~ 2 . 

The red line (4.7) intersects the q 2 axis at 

-1.8 • 10~ 5 , (4.8) 

in excellent agreement with the expected value q 2 = (3.23). The data implies 

cc Z c «* ex q" 2 <=► X ^o, q cx q~ 2+ " , (4.9) 

which determines the anomalous critical exponent rj as 

r] = 0. (4.10) 

Finally, we turn to the dynamical critical exponent z. According to (3.27), it can 
be extracted from the scaling of the pole in the dynamical susceptibility in the vicinity 
of the critical point (in the hydrodynamic limit) 



= y _1 



IXO* 



cx t~ 



(4.11) 



A typical behavior of the dynamical susceptibility in a thermodynamically stable phase, 
i.e., for p u = 0.0351 < p\ x (blue curve), and a thermodynamically unstable phase, i.e., 
for pn = 0.0353 > p\ x (red curve) is presented in Figure 6. We used q = 10~ 2 . Notice 
that in the stable phase, dynamical susceptibility has two separate poles 



= X m = 



ro=ro»,q=10 -2 ' 



(4.12) 



for ztt>* > — according to (3.25) these poles correspond to different relaxation time- 
scales r of the system at criticality. In the unstable phase two additional poles appear 
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Figure 7: (Colour online) The inverse of the dynamical susceptibility of M = 2* plasma 
in a stable phase (blue curve) and an unstable phase (red curve) at q = 1CT 2 for 
\ito\ < 1. 

for small values of itv. Figure 7 zooms in on the range of \ito\ <C 1 in the dynamical 
susceptibility. One of these poles is at negative values of ito*, corresponding to a 
negative relaxation time r. A detailed analysis show 7 that the negative relaxation time 
scales as r _1 = zro* oc — q. Clearly, negative relaxation times signal the instability in 
the system — this is precisely the instability due to the hydrodynamic (sound-channel) 
modes which propagate with c 2 s < once pn > p\ x , see Figure 1. Such instabilities 
are expected on general grounds: whenever a thermodynamic phase of a system has a 
negative specific heat, the hydrodynamic modes in the system are unstable [27] 8 . What 
we demonstrated here is that such instabilities also result in the negative relaxation 
time: i.e., instead of approaching the equilibrium a system is driven away from it. 
We now focus on poles in the dynamical susceptibility at ixv > 0: 



(4.13) 

ro={ro*, L (q) , ro*, fl (q)} 



where L and R are the indexes of the two positive poles, such that the corresponding 
relaxation rates 

(27rTr L (q))- 1 = ztD*, L (q) < (27rTr K (q))- 1 = ito*, R {q) . (4.14) 

We performed numerical analysis for different values q = {10~ 4 , 10~ 3 , 10~ 2 } — the 
poles itD* i {£,jj}(c|) have a well defined hydrodynamic limit q — > 0, which is obtained by 



7 Available from authors upon request. 

8 The reverse is not true: a thermodynamically stable system might still have instabilities [28]. 
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Figure 8: (Colour online) The relaxations rates (blue dots) itO^{L,R} ({bottom,top}) of 
M = 2* plasma in the vicinity of the critical point. The solid red lines are the quadratic 
fits to the data, and the dashed green line represents pn = p\ x . 

computing the susceptibilities ( at different values of pn ) strictly at q = 0: 

lim itt>^ {L ,R}(q) = itt>*,{L,R}(0) = zW*j L m ■ (4.15) 

The results of such analysis are presented in Figure 8. The blue dots are the relaxation 
rates of M = 2* plasma at criticality, the solid red lines are the best quadratic fits to the 
data. The top dots/curve corresponds to itu*^, and the bottom dots/curve corresponds 
to itn^L. Once again, the vertical green dashed line corresponds to p n = p\ x . The 
results of the analysis show that in the hydrodynamic limit, both the relaxation rates 
are finite 

(2ttT c r^})" 1 = ito %{L>R} ex (A Pll )° cx t° oc (2ttT c £)° ■ (4.16) 

Thus, 

r {L>/?} «f «(° z = 0. (4.17) 

In the previous section we pointed out that the critical behavior of M = 2* plasma 
should be identified with that of 'model A' according to dynamical critical phenomena 
classification in [25]. As such, the dynamical critical exponent z was predicted to be 
(3.16) 

z =2 + c-0 = 2, (4.18) 

prediction 
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which differs from the value we obtained (4.17). 
5 Conclusions 

In this paper, building up on the previous work [14-16], we presented a detailed analysis 
of the static and dynamic critical phenomena in strongly coupled M = 2* plasma. This 
model is a string theory derived example of gauge theory /gravity correspondence where 
one deforms M = 4 SYM by giving a mass m& to bosonic components and a mass to/ to 
fermionic components of M = 2 hypermultiplet. Generically, i.e., when to& ^ nif, such 
a deformation completely breaks the supersymmetry. At finite temperature M = 2* 
gauge theory plasma undergoes a second-order phase transition, provided m'j < m^. 
This continuous transition is characterized by a terminal temperature T c , which can 
be reached within isotropic and homogeneous equilibrium phases. At temperature 
T c the two phases continuously meet, with vanishing speed of sound. One of the 
phases is always perturbatively unstable — the instabilities reside in the sound-channel 
hydrodynamic modes which propagate with c 2 s < 0. Extending [15], we computed the 
static and the dynamical critical exponents of the transition, as approached from the 
perturbatively stable phase: 



As expected — since the gauge theory is a large- N model — the transition is of the 
mean-field theory type with vanishing anomalous static critical exponent r\. Under the 
static scaling hypothesis, only two of the critical exponents are independent — thus, 
the six static critical exponents in (5.1) must satisfy four algebraic constraints. Similar 
to analysis in [16], we find that only one of these constraints, i.e., the hyperscaling 
relation, is being violated 



Dynamical features of the transition in M = 2* plasma are quite interesting. First 
of all, the symmetries of the transition identify it as the one in the universality class 
of 'model A', according to classification of Hohenberg and Halperin [25]. The latter 
predicts the dynamical critical exponent as z pre( nction — 2, which contradicts direct 
computations (5.1). Both the stable and the unstable phases have multiple (two) 
relaxation times which remain finite at the critical point — hence the critical exponent 
z — 0. Once again, as in analysis in [15], even though the dynamical critical exponent 





2 - a ^ 3u. 



(5.2) 



20 



z 7^ 1, and thus there is an anisotropy between the time- and the space- coordinates 
scaling in the two-point (non-equilibrium) correlation functions, the dual gravitational 
geometry at criticality does not exhibit a Lifshitz-like scaling in the sense of [29]. In fact, 
if, as suggested by [16], different non-equilibrium correlation functions at criticality have 
different dynamical exponents z, it is not possible to 'by hand' embed the anisotropic 
scaling of the correlations functions into the symmetric of the dual geometry. The 
critical phenomena in [16] and the one considered here indicates that anisotropic scaling 
is rather an emergent phenomena. Second, the unstable phase of M = 2* plasma has 
an additional relaxation rate r~^ stable oc — \k\. A negative relaxation rate indicates that 
rather than approaching the equilibrium, a perturbed system is driven away from it 
- this is yet another reflection of the instability in the hydrodynamic sector of the 
theory, which is necessarily linked to a thermodynamic instability of the corresponding 
plasma phase [27]. 

In the future, it would be interesting to understand how the general classification of 
dynamical critical phenomena [25] should be enlarged to incorporate the universality 
class of M = 2* plasma. Probably the most pressing question is the understanding of 
the equilibrium phases in this universality class (and also the one of the M = 4 SYM 
plasma with an R-symmetry chemical potential) for T < T c . Such phases can not be 
homogeneous and isotropic. 
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